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Abstract 

The unitarity relations for isovector nucleon form factors in the two-pion approx- 
imation for the intermediate states are analytically extended in a model independent 
way from the real half-axis t > 4/j, 2 (/i is the pion mass) to the whole complex t- 
plane. As a result, the explicit integral representations for the nucleon form factors 
are constructed in terms of the pion form factor and imaginary parts of the t-channel 
p-wave 7r TV-scattering amplitudes at t < 0. New method based on the explicit integral 
representations is developed for numerical evaluation of the two-pion contribution to 
the isovector nucleon form factors. The results of different methods are compared and 
their efficiency is discussed. 
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1 Introduction 



Recent years new precise experimental data were collected for the electromagnetic 
nucleon form factors at low and intermediate momentum transfers at spacelike and timelike 
regions 0. The experimental progress stimulates development of the phenomenological 
models for the nucleon form factors. 

The unitarity relations for isovector nucleon form factors, discovered by Frazer and 
Fulco 0, appeared to be the first fruitful attempt to provide field-theoretical framework 
for description of the electromagnetic nucleon form factors at low momentum transfers. 
Being combined with the requirements of the analyticity, relativistic invariance, and crossing 
symmetry of the tt ^-scattering amplitudes, the unitarity relations allow to evaluate the two- 
pion contribution to the absorptive pars of the nucleon form factors at the cut (4/i 2 , +oo) 
using experimental data on the 7T7r-and 7riV-scattering amplitudes only. Similar relations 
based on the unitarity are used in the analysis of the electromagnetic pion form factor M. 
The unitarity relations explain the origin of the vector-meson dominance (VMD) in the pion 
and nucleon form factors. The VMD models are in wide use in the modern calculations [HI 
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The differential approximation scheme for the analytical continuation of the ttN- 
scattering amplitudes from the s-channel (ttN — > ttN) to the t-channel (tttt — > NN), sug- 
gested by Efremov, Meshcheryakov and Shirkov [|l^], was used by Lendel et al. |L3| to 



evaluate the nucleon form factors from the unitarity relations (for a review see Ref. ||14|| ). 



Hhler and Pietarinen JT5| performed the alternative analysis using the prescription by 
Frazer and Fulco || for calculation of the ^-channel 7riV-scattering amplitudes at t < and 
the discrepancy function method to extrapolate the amplitudes to the region t > 4/i 2 (for a 
review see Ref. [pf]). 

The agreement of these two methods is surprisingly pure. The estimates of the two- 
pion contribution to the isovector nucleon form factor F 2v (t = 0) are in reasonable agree- 
ment, whereas the estimates of the form factor F lv (t = 0) differ by a factor of two. There 
were no attempts to analyze the origin of this discrepancy. This might be a reason for the 
unitarity relations are simply ignored during the last two decades in the most calculations 
of the nucleon form factors. 

The apparent disagreement of these two methods indicates presence of uncertainties 
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whose origin is not yet understood clearly. The purpose of this work is to discuss possible 
sources of the errors and develop new method free from ambiguities inherent in the earlier 
calculations of isovector nucleon form factors. 



The unitarity relations by Frazer and Fulco are derived in the two-pion approximation 
for the intermediate states between the photon and the iViV-pair (see Fig.l). Integral rep- 
resentations for the nucleon form factors, equivalent to the FF unitarity ||, are constructed 
in Refs. |L7], [R| using the methods of Refs. [[|, |12|-[p4| and Refs. || [15], [UJ and assuming 



the validity of the Frazer-Fulco-Gounaris-Sakurai (FFGS) model for the pion form factor 
fH, |l9fl . These representations extend the FF unitarity from the cut (4/i 2 , +oo) to the com- 
plex t-plane. For the kaon form factor, the representation of such a kind is constructed in 



Ref. 1 20]. In this paper, we extend the FF unitarity to the complex t-plane in a model 
independent way and give numerical estimates for the nucleon form factors, based on these 
represent at ions . 

The outline of the paper is as follows. In the next Sect., we construct analytical repre- 
sentations for the nucleon form factors valid in the complex t-plane, completely equivalent 
to the FF unitarity. In Sect. 3, we develop new method for numerical evaluation of isovector 
nucleon form factors, based on these representations. In Sect. 4, we compare our results with 
results of the previous works and discuss efficiency of different methods for evaluation of the 
two-pion contribution to the nucleon form factors from the unitarity. 



2 Analytical Extension of the Unitarity Relations 

The isovector nucleon form factors are defined by isovector part of the electromagnetic 
nucleon current 

3n{p',s',P,s)v = u{p' ,s')[F lv (t)^^ + — F 2v (t)ia IJ , l/ q l/ }T 3 u(p, s) (2.1) 

where t = (p' — p) 2 , m the nucleon mass, and r 3 the third isotopic Pauli matrix. The form 
factors are normalized by Fi v (0) = 1/2 and F 2v (0) = (/i p — fi n — l)/2 with \x v and \i n being 
the proton and neutron magnetic moments (in n.m.). The proton charge is set equal to 
unity. 

The unitarity relations for nucleon form factors have the form M 

ImF iv (t) = ~F*(t)Ti(t) (2.2) 
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where k = yt/A — /i 2 , /i is the pion mass, and 

m^^i-fi+A^f 1 -}, 

pi 4V2m 

with p 2 = m 2 — t/4. The are the t-channel p-wave 7riV-scattering amplitudes. The 

unitarity relations are shown graphically in Fig. 1. 

In the approximation of neglecting all but the two-pion intermediate state, the values 
Ti(t) obey the dispersion relations 

/ \ / 7—1 / \ 1 f a ImTAt') dt' . , 

Ti(t) F«(t) = - / - l \ } 2.4 
7T J-00 V — t r w {t') 

with 
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a = 4/i 2 - ^. (2.5) 



Eqs.(2.2) and (2.4) are exactly valid below the four-pion threshold at the interval 



4/i 2 < t < 16/i 2 . The branching ratio |^TJ B(p — > Ait) < 0.002 is small, so the two-pion 
approximation for the unitarity relations is valid with high accuracy at the p-meson peak 
and even at higher energies. 

It is useful to define a function 

K * {t) = nL hT* V* — m)l (2 - 6) 

with t\ > 0. This function is analytical in the complex t-plane with a cut (4/i 2 , +00). The 
factor ti/(t\ + 1') is introduced to ensure better convergence of the integral. With the help 
of the function (2.6) the unitarity relations (2.2) take the form 

ImF m (t) = -^^ImK^Ti^/F^t). (2.7) 

Using Eq.(2.4), the once subtracted dispersion integral for the nucleon form factors can 
easily be evaluated to give 

t\ 7T J-00 t — t -Tn{t ) 

For computation of the nucleon form factors, it is sufficient to know the imaginary parts 
ImTiit) at t < a. The integral in Eq.(2.8), however, converges very slowly. It is determined 
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by low- and high-energy parts of the values ImTiit). The high-energy part —oo<t<b 
is unknown. It has been shown that there are no objections against the truncated 
partial wave expansion in the region - 40fi 2 < t < 0, so we set b = — 40/i 2 . The values 
can be decomposed to the Born and rescattering parts: Tj(t) = T iB (t) + T iR {t). The real 
and imaginary parts of the values are know on the interval (6,0) from the phase- 

shift analysis of the 7riV-scattering amplitudes. The problem consists in an approximate 
evaluation of the high-energy part with the help of the information on the low-energy real 
pars ReTi(t). 

The once subtracted dispersion representation for the values Ti(t) has the form 

The value k is assumed to be real, k < a. The analog of Eq.(2.8) reads 

F iv (t) = F iv (-h) + ' 



t 



r„ r<(*c) , l n f a ImTiWit-KtKM-V - k)^) dt' 

[K6 f 7T (k) + nJ-oo (f-K) (t-f) F n (t'y- 1 Uj 

In comparison with Eq.(2.8), the integral at high t is suppressed by the additional power of 
t. Eqs.(2.8) and (2.10) constitute the main result of this section. 



These equations can be simplified when the FFGS model 0, [L9| for the pion form 
factor is applied. In this model, the pion form factor has the form 

where 

D(t) = A + Bt + k 2 h{t) (2.12) 

is the D-function. It determines the phase shift 8\i(t) in the I = J = 1 7T7r-partial wave. The 
exponentially large parameter t oc /x 2 exp( ) is a unique root of equation D(— to) = 0. 
Here m p and T p are the p-meson mass and width, k p = ^Jm^/A — /i 2 . The parameters A 
and B are fixed by requirements 

Mm 2 )=7r/2, 
6' n {m*) = l/(m p T p ). 

In this way one gets 

A = m 2 p + n 2 h p + m 2 p k 2 p h' p , 
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B = -1 - \h p - k 2 p h' p . (2.13) 



The function h(t) is determined at t < by expression 



K , U \ 1 /V-*, / W-* + v 7 ^ 



m 2 p T p w 2tt V -* VV-^-v^ 

The analytical continuation of this function to the complex t-plane has a cut at t > 4/i 2 . At 
the interval < t < 4fi 2 the function h(t) is analytical and real. At t > A/i 2 

R e h(t) = t- 4//2 fe(t- 4// 2 ), 

A; 3 fc 
p Imh(t + iO) = ■=. (2.15) 



In Eqs.(2.7) /i p = Reh(m 2 ) and ft,' = Reh'(m 2 ). The D-function at the origin takes the 
form 

D(0) = m 2 p + ^h p + m 2 p k 2 p h' p - (2.16) 

In the narrow width limit, A = m 2 , B = —1, D(0) = m 2 , t = oo, and so the F n (t) takes a 
monopole form. 

The following relation takes place above the two-pion threshold 

4l^)| 2 = D(0) % + iO). (2.17) 

vt m p r p t 

The dispersion integral (2.6) for t\ = to can then be evaluated to give 

K n (t)=D(0)^-F n (t). (2.18) 
m 2 p T p 



Eq.(2.10) becomes p 



, TJk) 1 /"« ImTAt'){t- K)/FJt') - it' - K)/FJt) ,,, 

(Re^-{ + -P / 11 n / \ 7rV J dt'). 2.19 

v F v {k) tt 7-oo (t'-K) (t'-t) ' K J 

The imaginary part of the ratio F iv (t) j 'F^(t) has a simple representation: 

Im{(F iv (t) - F lv (-t ))/F w (t)}/^= = -- f a I -^p-dt'. (2.20) 

y/t VT J-oo t — t 
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3 Numerical Evaluation of the Nucleon Form Factors 



In this Sect., we give numerical estimates for the nucleon form factors using the FFGS 
model for the pion form factor. 

In the real part of the ratio F iv (t)/F n (t), convergence of the integral can be improved 
by taking a weighted sum 



Re{(F lv (t) - F tv (-t ))/F n (t)} = -D(0) 



kp t + t 
m 2 V p t 



VV (P M K i) i l P f a ImT ^ - K i)l F «W) - (g - ^)Re{l/F w (t)} 
with the coefficients 

N 

E% = 1 - (3-2) 

The subtraction points Kj are assumed to belong to the interval b < kj < 0, so the 
values Re yj^j entering Eq.(3.1) are known. The coefficients c^- can be chosen such as to 
suppress the high-energy part — oo < t < b of the integral. 

The integral in the right hand side of Eq.(3.1) can be rewritten in the form 

N r< I., \ 1 ,h rn. N 



3=1 

where 



E + -( / +P / ) E c l3 ImT t {t')^ f , Kj )dt? (3.3) 

~[ ^tt(«j) TT J-oo Jb ■ \ 



, v 1 (t - K 3 )/F^t>) - (f - K 3 )Re{l/F^t)} 

The first and third terms are known. We treat square of the second term as an error to be 
minimized. It can be evaluated as 

I rb N 1 rb N 

(- / E CijlmTiit')^, t', K^dt'f « (- / E CijlmFiB^it, t', n^dt'f < 

7T J-oo - =1 7T J-oo - =1 

1 /•& r-b N b 3 

- {ImT iB {t')fdt' (Y. c *Mt,t',K J )) 2 dt' = el/(D(0)^-) 2 (3.5) 

7T J— oo J — oo =1 ^p-L p 
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The second source of the error is finite precision of the rescattering part of the values 
ReT iR (t). The corresponding contribution to the error of the expression (3.3) reads 



P N r (ic \ 

4 = P(o)^) 2 E4(^^tt) 2a ? ( 3 - 6 ) 

where Aj = 0.05 are relative errors for the ReTmit). When N goes to infinity, the values 
Aj are no longer statistically independent, in which case expression (3.6) overestimates the 
error. In this work, we restrict ourselves by relatively small values N < 12. The errors of 
similar nature coming from the finite precision of the ImTm(t) are neglected. 

The conventional minimum of the total error e\ ot = e\ + e| with respect to the coef- 
ficients Cij under the constraint (3.2) imposed can be found by introducing the Lagrange 
multiplier A. The coefficients Cy and A are solutions of a system of the N + l linear equations 
obtained by taking derivatives with respect to the c^j and A of the function 

TV 

el + el-X^c^-l). (3.7) 

The coefficients c^- and A are t-dependent. 

The real and imaginary parts of the rescattering amplitudes Vm on the interval (b, 0) 



are calculated in Ref. |H| and more recent results can be found in Ref. fl6 |. 



We test first convergence of the method with increasing the number N of the weight 
coefficients c^-. In Figs. 2(a,b) the values of the form factors at zero momentum transfer 
are shown, together with the errors e to t, as functions of the number N. The subtractions 
are made at the points Kj = b(j — 1/2) /N. For N > 6 the errors decrease slowly. The 
predictions are more stable for the F 2v (0). It is seen that variations of the Fi v (0) are within 
the error bars. The variations of the ^^(O) are smaller than the theoretical errors. The 
coefficients have irregular behavior with increasing the N. The results shown in Figs. 
2(a,b), however, agree with each other for different values of the N. This consistency check 
is equivalent to the consistency check for the analytical continuation of the amplitudes Tj 
to the region t > Ajj 2 with the use of the discrepancy method. 

The convergence of the integral in Eq.(2.20) for imaginary part of the ratio F iv {t) / F n {t) 
can be improved by adding to the right hand side a weighted sum of expressions 

TAkA 1 f a ImTAt') dt' , . 

Re *\ J { - -P / tK ' - . = 0. 3.8 
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The integrals in Eqs.(2.20) and (3.8) have similar structure. We chose the weights dij to 
suppress the high-energy part of the integral (2.20). 

The right hand side of Eq.(2.20) takes the form 
N T(k) 1 r b r a N 1 1 1 

~ * ( L +P L wM^ir^m + —? d< < 3 - 9 > 

To ensure faster convergence of the integral, we require 

N 

J2dij = 0. (3.10) 

3=1 

The integral can be regularized afterwards by a substitution l/(t' — Kj) — > l/(t' — Kj) — 

l/^-Ko). 

The error coming from the unknown high-energy part of the integral can be evaluated 

to be 

« = h Ij^iW fjE ^ (f _ %) V_ Ko) ^ + (3-n) 

The error coming from the finite precision of the rescattering parts ReTm(t) reads 

^E4(^#f)^ 2 (3-12) 

j=l r ^\ l% j) 

The conventional minimum of the total error S^ ot = 8f + 8% can be found with the help 
of the Lagrange multiplier 7 corresponding to the constraint (3.10). The coefficients and 
7 are solutions of a system of the iV + 1 linear equations obtained by setting equal to zero 
derivatives with respect to the d^ and 7 of the function 

N 

81 + 81-^^3- (3-13) 

j'=i 

In Figs. 2(c,d) the values Im{(Fi V (t) — Fi V {— to)) / F n (t)} at the p-meson peak t = m 2 p , 
together with the errors, are shown as functions of the number N. The results seem to 
converge with increasing the N. The variations are within the error bars. The predictions 
are more stable for the F 2v {mfy than for the F lv (m 2 p ). The variations of the F 2v (m 2 p ) are 
noticeably smaller as compared to the errors. 
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The theoretical errors are defined by replacing — > IYb(£) in the integral over the 
region (—00,6). These estimates are valid up to unknown coefficients. These coefficients 
can apparently be fixed by requiring that the theoretical errors be equal to the statistical 
ones. In case of the form factor F lv (t), such a coefficient is close to unity, whereas in case 
of the form factor F2 V (t) such a coefficient should be 5-7 times smaller than the unity. 

Shown in Figs. 3(a-d) are real and imaginary parts of the ratios between the nucleon 
form factors and the pion form factor versus t. The errors etot for the real parts and Stot for 
the imaginary parts are shown. The visible t-dependence of the errors Stot i n Figs. 3(c,d) is 
mainly due to the logarithmic scale of the plots. 



4 Discussions 



In the differential approximation scheme by Efremov, Meshcheryakov and Shirkov [12 
the t-channel p-wave projections of the 7riV-scattering amplitudes are calculated with an ap- 
proximation of the p-wave amplitudes by linear superpositions of the forward and backward 
amplitudes and their higher derivatives. The analytical continuation to the unphysical re- 
gion t > Afi 2 is performed with the help of dispersion relations for the forward and backward 
amplitudes. On the basis of the unitarity relations by Frazer and Fulco H , the nucleon form 
factors are analyzed to the lowest order of the differential approximation by Lendel et al. 
II (see also Ref. 0). 



In the original work by Frazer and Fulco, the following program was suggested for 
evaluation of the nucleon form factors from the unitarity relations: The invariant s-channel 
7riV-scattering amplitudes at fixed t < are continued analytically from the physical region 
to the unphysical region between the s- and w-channels. Afterwards, the t-channel partial 
wave projections of these amplitudes can be computed. The complex p- waves f±(t) are 
known in the region - 26 fi 2 < t < in which the partial wave expansion converges. It has 
been shown [JZ2] that there are no objections against the truncated expansion also in the 
region - 40/i 2 < t < 0. The analytical continuation of the complex ]>waves f±(t) to the region 
t > Afi 2 can be made with the use of the extended unitarity relation, Eq.(2.4), according 
to which the ratios f±(t)/F n (t) (equivalently Y i{t) / F w {t)) are analytical functions in the 
complex t-plane with the left branch cut (—00, a). In practice, the discrepancy function 
method is suitable in a restricted region of the complex t-plane for solving this task, since 
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the high-energy behavior of the scattering amplitudes is unknown. Given that the analytical 
continuation to the region t > 4/z 2 is performed, one can compute the dispersion integral for 
the nucleon form factors. As t increases, the analytical continuation becomes less precise, 
so the upper limit A 2 « 1 GeV 2 should be introduced to exclude the high-energy part of 
the dispersion integral. This program was performed, within the framework of the FFGS 



model, by Hoehler and Pietarinen |15j for the detailed analysis of the isovector nucleon form 
factors (see also Ref. (1161). 



In the latter approach, one should compute three integrals numerically: The first one 
comes from the computation of the t-channel partial wave projections of the ttN amplitudes 
at t < 0. The second one comes from the dispersion integral for analytical continuation of 
the p-wave amplitudes f±(t) to the region t > 4/r 2 . The third one comes from the dispersion 
integral for the nucleon form factors. 

We demonstrated that the analytical evaluation of the last integral is also possible. In 
doing so, we arrived to the explicit one-dimensional integral representation (2.19) for the 
nucleon form factors in terms of the pion form factor and imaginary parts of the amplitudes 
Fi(t) at t < a. There is no therefore need to perform analytical continuation of the ampli- 
tudes Ti(t) to the region t > Ap 2 . The representation (2.19) is constructed for the FFGS 
model. The model independent integral representations are given by Eqs.(2.8) and (2.10). 

In Eq.(2.8), the imaginary part of the amplitudes Tj(t) occurs only. Due to the ex- 
tended unitarity (2.4), the low-energy real parts ReTi(t) are expressible through the high- 
energy imaginary parts ImTi(t). Therefore, the real parts contain information about high- 
energy behavior of the amplitudes. This information can be used for evaluation of the 
high-energy part of the integrals. In Sect. 3, we described the corresponding numerical algo- 
rithm. 



To get an idea about accuracy of the method |T2j - [|14|], one needs to calculate the 
next order term in the differential expansion. This work is not done yet. It is clear, however, 
that the Born amplitude should be treated beyond the differential approximation method, 
since the logarithmic singularity at t — a provides pure convergence of the series expansion. 



The errors of the analysis [15j plotted on Fig.l are ours. They are calculated using the 
claimed accuracy ±15% in Ref. ]T5| for the amplitudes f±(t) at the p- meson peak. There 
exist additional errors not plotted on Fig.2. 
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(i) The existence of the cut-off parameter A 2 m 1 GeV 2 implies that a contribution of 
the p-meson tail t > A 2 to the nucleon form factors is disregarded. The results [15| are A 
-dependent. In our approach, the cut-off parameter is safely sent to infinity (the integrals 
converge). 

(ii) In the analysis fL5 |, the contribution from the interval -120/i 2 < t < b, in which the 



truncated partial wave expansion of the amplitudes diverges, is included in the dispersion 
integral. The partial wave expansion gives at the interval -120/i 2 < t < b substantial part of 
the discontinuity ImTi(t). The corresponding errors cannot be controlled, since the values 
ImTm(i) calculated from the partial wave expansion exhibit at t < b strong tendency to 
increase with decreasing t. 

The contribution from the interval (—00, —120/i 2 ) in Ref. [pl| is interpolated by a 
polynomial a + bt. From the unitarity it follows that Ti(t) = 0(1 /t) and T 2 (t) = 0(l/t 3 ^ 2 ) 
as t I— > 00. Taking into account that in the FFGS model F w (t) = 0(l/log(t)), we conclude 
that the dispersion integral in Eq.(2.4) can be evaluated, respectively, as 0(log 2 (t) /t) and 
0(1/ 1) as t 1— > 00. The contribution from the interval (—00, —120/i 2 ) therefore vanishes when 
t increases. The polynomial representation a + bt does not satisfy this requirement. The 
simple interpolations a/(b + t) and log((a + t)/(b + t)) could provide the correct asymptotic 
form. 



The distinction between the methods |L2] - |L6| and ours makes comparison of the 
results not so straightforward. 

The comparison of the results makes sense for zero value of the subtraction constant 
Fi V (—to) = 0. The quark counting rules imply that the total nucleon form factors, containing 
all multipion and AW contributions to the nucleon spectral functions, vanish as t goes to 
infinity. In the two-pion approximation, the value F iv (—to) should not, however, be strictly 
equal to zero. 

In Fig. 2(a-d), we show our results at two points t = and t = m 2 for different 
numbers of the weight coefficients. The errors are theoretical ones, etot and (k 3 /\/i)5, 



tot- 



The results of Ref. [15] are also displayed. 



The most typical results are summarized in Table 1. Ref. |15j gives for the amplitudes 
f±(t) at the p -meson peak an error of ±15%, resulting in turn in a small error of the form 
factor Fi v (0), which is 4-5 times smaller as compared to our predictions. As we already 



12 



discussed, some of the uncertainties in the analysis of Ref. can hardly be controlled, so 



the errors of Ref. |0|] are greater as compared to those shown on Figs. 2 and 3 and in Table 
1. Our method for large N gives systematically lower (negative) values for the difference 
Fi v (0) — F lv (—to). The error e to t is, however, large. The unitarity relations, therefore, do 
not allow to give accurate predictions for the form factor F\ v (t). The value -F2„(0) in our 
approach is well defined. The estimates of Refs. [jn| and |15| are in reasonable agreement 
with ours. The accuracy of our predictions for the -F 2? ,(t) is better as compared to the earlier 
ones. 



5 Concluding Remarks 

The unitarity relations by Frazer and Fulco are written in the two-pion approximation 
for imaginary parts of the isovector nucleon form factors at the branch cut (4/i 2 , +oo). 
We showed that the unitarity relations for the nucleon form factors can be solved to give 
analytical one-dimensional integral representations for the nucleon form factors valid in the 
whole complex t-plane (Eqs.(2.13) and (2.15)). The nucleon form factors are expressed in 
terms of the pion form factor and the t-channel p-wave amplitudes Tj(t) at t < a. 

The representation (2.15) is used to develop new method for numerical evaluation of 
the nucleon form factors from the unitarity relations. The method is described in Sect. 3. It 
has some advantages as compared to the earlier methods: (i) there is no need to perform 
analytical continuation of the amplitudes Tj(t) from the region t < to the region t > 4/i 2 , 
(ii) no cut off parameter A to be introduced in the dispersion integrals for the nucleon form 
factors. 

The numerical estimates of the nucleon form factors from the unitarity relations in 
the two-pion approximation are at present in the rough qualitative agreement for different 
methods. These estimates are considered as efficient and reliable if results of the different 
methods will be brought into the correspondence by evaluation of the high-order approxi- 
mations and careful analysis of possible sources of errors. 
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FIGURE CAPTIONS 



Fig. 1. Pictorial representation of the unitarity relations for isovector nucleon form 
factors in the two-pion approximation for the intermediate states between the photon and 
the nucleon-antinucleon pair. The imaginary parts of the nucleon form factors are expressible 
in terms of a product of the pion form factor and the t-channel 7riV-scattering amplitudes. 

Fig. 2. (a,b) The values F iv (0) — F iv (—t ) (left scale) versus the number N of the 
weight coefficients fixed by a requirement of minimum of the total errors e to t originating 
from the unknown high-energy part of the integrals and finite precision of the real parts of the 
t-channel p-wave scattering amplitudes. (c,d) The ratios Im{(F iv (m^) — F iv (—t ))/ F n (mfy} 
(right scale) at the p-meson peak versus the number N of the weight coefficients fixed by 
a requirement of minimum of the total errors 5tot originating from the unknown high-energy 
part of the integral in Eq.(3.9) and finite precision of the real part of the t-channel p-wave 
scattering amplitudes. The results of Ref. (dashed regions) are plotted for comparison. 

Fig. 3. (a,b). Real part of the ratios between the nucleon form factors and the pion 
form factor, Re{(F iv (t) — F iv (— t ))/F n (t)}, versus t with errors e tot for N = 6 (left scale). 
(a,b) Imaginary part of the ratios between the nucleon form factors and the pion form factor, 
Im{(F iv (t) — Fiv(— to)) / F^(t)} / versus t with errors 8 tot for N = 6 (right scale). The 
results of Ref. |l5j (dashed regions) are plotted for comparison. 
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TABLE 1. Comparison of predictions of different methods at two points t = and 
mi. In Ref. ||13|| , the estimates are obtained in the first order of the differential 



approximation. The errors (presently unknown) can be estimated by computing the next 
order of the differential approximation. The errors of Ref. |T5[ are evaluated from the 
claimed accuracy ±15% for the amplitudes f±(mi). Two sets of values are given on the 
basis of our calculations: shown in the line (I) are values corresponding to the maximum 
number of the weight coefficients N = 12. In the line (II) predictions obtained by statistical 
averaging the data from N = 3 to 12 are given. In the line (I), the errors are theoretical 
ones, e tot and (k 3 /Vi)5 tot (see the text), while in the lower line the errors are statistical ones. 
The reduction of the errors for the case of the form factor F2 V (t) occurs, since the variation 
of the results for different N in this case is significantly smaller than the theoretical errors, 
e to t and (k 3 /y/i)5 t 



Hot- 



Ref. (F iv (0) - F iv (-t o ))/F n (0) Im{(F iv (m 2 p ) - F iv {-t )) / F^rn^} 

i=l 2=2 i = 1 i = 2 

|3| 1.15 3.09 0.05 0.09 

PI 0.52 ±0.08 3.15 ±0.57 0.20 ± 0.09 1.00 ± 0.21 

This work (I) -0.36 ± 0.35 2.73 ± 0.50 0.39 ± 0.07 1.08 ± 0.07 

This work (II) -0.03 ± 0.37 2.78 ±0.12 0.34 ±0.05 1.07 ±0.01 
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